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in 

Q ' Abstract 

We present a summary of work done on dense hadronic matter, based on the Skyrme model, 
which provides a unified approach to high density, valid in the large Nc limit. In our picture, 
■ dense hadronic matter is described by the classical soliton configuration with minimum energy 

for the given baryon number density. By incorporating the meson fluctuations on such ground 
Qh! state we obtain an effective Lagrangian for meson dynamics in a dense medium. Our starting 

Q_i' point has been the Skyrme model defined in terms of pions, thereafter we have extended and 

1^ ■ improved the model by incorporating other degrees of freedom such as dilaton, kaons and vector 

mesons. 



X 



1 Introduction 

At high temperature and/or density, hadrons are expected to possess properties that are very 
different from those at normal conditions. Understanding the properties of hadrons in such ex- 
treme conditions is currently an important issue not only in nuclear and particle physics but also 
in many other related fields such as astrophysics. Data from the high energy heavy ion colliders, 
astronomical observations on compact stars and some theoretical developments have shown that 
the phase diagram of hadronic matter is far richer and more interesting than initially expected. 
Lattice QCD calculations have been carried out successfully at high temperature, however sim- 
ilar calculations at high density have not yet been possible^. Theoretical developments have 
unveiled such interesting QCD phases as color superconductivity Hj. Moreover effective theories 
can be derived for these extreme conditions, using macroscopic degrees of freedom, by matching 
them to QCD at a scale close to the chiral scale ~ ^T^fn ~ 1 GeV[H]. 

We have followed a different path to dense matter studies by using as our starting point a 
model Lagrangian, in the spirit of Skyrme, which describes hadronic matter and meson dynamics 
respecting the symmetries of QCD. The parameters of the model are fixed by meson dynamics 
at zero baryon number density. A la Skyrme|5|, baryons arise from a soliton solution, the 
skyrmion, with the topological winding number describing the baryon number. In our scheme 
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dense matter is approximated by a system of skyrmions with a given baryon number density 
whose ground state arises as a crystal configuration [Sj 1^1 • Starting from this ground state our 
approach provides insight on the intrinsic in-medium dependence of meson dynamics. We have 
studied (i) the in-medium properties of the mesons and (ii) the role of the other degrees of 
freedom besides pions in the description of matter as it becomes denser [71 151 1^ ITUl ITTl [T^ . 



2 Model Lagrangians 

The original Skyrme model Lagrangianjl] reads 

= -llTTiL.m + ^Tr[L^, L^f + ^Tt{U + C/t _ 2), (1) 

where = U^d^U and U = exp{iT ■ vr) G SU{2) is a nonlinear realization of the pion fields, 
/tt the decay constant and the pion mass . The second term with e, the Skyrme parameter, 
was introduced to stabilize the soliton solution. 

The dilaton field x was incorporated in the model to make it consistent with the scale 
anomaly of QCD 13 . The Lagrangian then becomes 

= -f (t) Tr(L,L^) + ^Tr[L„L.]2 
+\d^xd''x-V{x), 

Note the different powers of {x/ fx) ™ front of each term. The last line is the Lagrangian for the 
free dilaton field, where V{x) = (m^/^/4)((x//x)^(ln(x//x) ~ i) ~ i)' "^x dilaton mass 

and its decay constant. 

The vector mesons, p and uj, can be incorporated into the Lagrangian as dynamical gauge 
bosons of a hidden local gauge symmetry which requires the doubling of the degrees of freedom 
as U = ^[^R. This La grangian reads (14j 

L.xP^ = -f(f) Tr(L^^^) + ^ (^-| j Tr(C/ + C/t-2) 

_ila (j^j Tr[CLd,C[ + ^Rd.^l + i(5/2)(f • p, + co,)f (3) 

with p^^ = df^piy — d^p^ + gpfj_ x p^, uj^^y = d^uj,^ — dyOJ^j,, and where i?^ is the topological baryon 
number current. The quartic Skyrme term of is not present, because its stabilizing role is 
played here by the vector mesons. 



3 In-Medium Pion Dynamics 

In order to explain our basic strategy, we start with the Skyrme model Table 1 serves to 
generalize the discussion to other models. The vacuum solution of in the i? = sector is 
U = \. Fluctuations on top of this vacuum describe pion dynamics. The Lagrangian supports 
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Figure 1: This figure shows the energy per skyrmion and the average value of cj = (^Tr(C/)) 
as a function of the FCC crystal parameter L (in units of (c/tt)"^). The figures in the small 
boxes correspond to samples of the baryon number distribution on the xy-plane {z = 0). The 
corresponding FCC boxes are drawn by a square. 

solitons with nontrivial topological structures. The B = 1 soliton solution with lowest energy 
is a "hedgehog" , U^^^ = exp(ir • fF{r)) with the boundary condition for the profile function 
F{r), -F(O) = TT and F{oo) = 0. The two-skyrmion system has lowest energy when one on 
the skyrmions is rotated relatively to the other in isospin space, by an angle tt, about an axis 
perpendicular to the line joining their centers. Therefore, skyrmion matter has the lowest energy 
in an FCC (face centered cubic) single skyrmion crystal structure, where the skyrmions at the 
nearest site are relatively oriented in such a low energy configuration and U (r) has the symmetry 
structure described in Table 1. 

For a given baryon number density, the lowest energy configuration can be found numerically 
by varying the field values at points of a discrete mesh^ or by adjusting the coefficients of the 
Fourier series expansion for the fields with specific symmetries 0. 

We show in Figure 1 the energy per baryon E/B of skyrmion matter and cr, the average 
value of ^TrC/ = (pQ over space, as a function of the FCC single skyrmion crystal size parameter 
L. Dimensionless units are used for E/B (in units of QTr'^f-j^/e) and L(in units of (e/7r)~^). For 
massless pions, for a value of L ~ 3.8 the slope of E/B changes indicating that the system 
undergoes a first order phase transition. If we look into the actual baryon number distributions 
(see figures in the small boxes), we see that the system changes from the FCC single skyrmion 
crystal to a CC half-skyrmion crystal. In the former phase, a well-localized single skyrmion is 
located at each FCC lattice site where (pQ = —1. In the latter phase, one half of the baryon 
number carried by the single skyrmion is concentrated at each FCC site while the other is 
concentrated on the links where (j)^ = +1. Both "half-skyrmions" centered at the points where 

= ±1 can hardly be distinguished and therefore they form a CC crystal. This phase transition 
can be seen more apparently in the quantity (^TrC/) (see the inset figure). In the literature jS], 
the vanishing of a is often interpreted as the restoration of the chiral symmetry. 
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If we turn on the pion mass (results shown in the figure by open circles), we see no sudden 
change in the slope of E/B but at sufficiently high baryon number density skyrmion matter 
behaves as in an approximate half-skyrmion phase. 

The phase to the left of the minimum, referred in our work as "homogeneous" , is described 
by a crystal configuration. The phase to the right of the minimum, which we called "inhomoge- 
neous", because the pressure P = dE/dV is negative and therefore skyrmion matter is unstable 
against condensation of the skyrmions into dense lumps leaving large volumes of space empty. 

We incorporate the pion fluctuations through the Ansatz 



U{r,t) = VU^UM{f)VU^, (4) 

where Um{^ describes the background baryonic matter and U-,^ = exp{iTaTi'a/ fw) describes the 
fluctuating pions. Substitution of this Ansatz into leads to 



(5) 



where we have expanded up to the second order in the fluctuating fields. Eq.@ describes 
the dynamics of the pion in the dense medium. Gabi^, ^(^) ^.re the interaction 

potentials appearing due to background matter. 

As a first approximation let's average the potentials over space. Due to the reflection sym- 
metries of the background matter, {Vp = (A^ = 0, {Gab) = (1 - <Aa0fe ) = Z'^^ab and {(Pq) = a. 
Thus, the Lagrangian becomes 

L^^o = IZ^d^n'^d^'^a + ^mla^a^a. (6) 

The factor in front of the kinetic term can be absorbed into the renormalization of the pion 
fields as 

TTa — > vr' = ZT^TTa- (7) 

In terms of this newly defined fields vr^, the Lagrangian can be rewritten as 

L^^' = hd,<d>^< + '^ml{a/Z^^)n'aK- (8) 

This implies that the pions effective mass in the medium is given by 

m* = m^^/Z^. (9) 

Furthermore, we may reinterpret the wave function renormalization factor Z^^ as the ratio be- 
tween the in-medium decay constant and the free space one, 

f:/fn = Z^. (10) 

The results on these in-medium quantities @ and H10|l are presented in Figure 2 as a function 
of the baryon number density. The scales are strongly dependent on the the parameters and 
e, e.g. the density units in p could change considerably, and therefore the shown numbers 
should not be taken as definitive. However, the qualitative behavior will remain unchanged. In 
particular, the dependence on the baryon number density showing approximate scaling with the 
Skyrme parameter e as ~ is a solid statement. 

In Figure 2, we show that at low densities m*^/m-,^ ~ 1, while at higher density the ratio 
decreases down to zero. The density dependence does not appear linear in p, since the classical 
background leads automatically to higher powers in the density dependence. At higher densities, 
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Figure 2: We show the effective pion mass and decay constant in the dense medium as a function 
of the baryon number density. The dependence on the baryon number density approximately 
scales with the Skyrme parameter e as ~ e~^. The right figure illustrates a pseudo-gap scenario. 

powers greater than one in p come to play important roles. As the density increases the decay 
constant /* decreases, which can be interpreted as a partial restoration of the chiral symmetry 
in the medium. However, the ratio f^/fn only decreases up to ~ 0.65 remaining constant 
thereafter. This non vanishing of the pion decay constant, despite the vanishing of a, indicates 
that we may be describing a phase which is not in the standard Wigner- Weyl symmetry. It may 
be a pseudo-gap phase where the gap is non-zero though chiral symmetry is restored, resembling 
what might be happening in the normal phase of high Tc super conductivity TBI. The pseudo-gap 
phase is schematically illustrated in Figure 2. 

The pseudo-gap phase may be an artifact of the model with only pions. If we have only pion 
degrees of freedom which are realized non-linearly through the phase of U, they must live on 
the chiral circle and the pseudo-gap phenomenon might appear. 

However if we introduce a dilaton field X) as in we see that it may shrink the chiral 
circle to a point and the conventional Wigner- Weyl symmetry appears 0. To illustrate this we 
take the x field constant, x/fx ~ ^- baryon free space the vacuum of the dilaton field is 
the minimum point of V{x) and is X = 1. However, for the dense baryonic matter, the ground 
state should be the minimum of E/B, the energy per baryon expressed as 



Here, E2, -E4 and Em are the contributions from the skyrmion field configurations of the given 
baryon number density through the kinetic term, Skyrme term and mass term of the Skyrme 
Lagrangian, respectively. Couplings to dilaton matter terms contributes with other effective 
potentials which will modify the value of X at which E/B has a minimum further. 

Shown in Figure 3 is E/B as function of X for various FCC length parameters, where we 
have substituted for E2/B, E4/B and E^/B the values that led to Figure 1. At low density 
(large L), the minimum of the effective potential E/B{X) is shifted slightly away from X = 1. 
As the density increases, E/B(X) starts developing another minimum at X = which was an 
unstable extremum in free space. At L ~ 1 fm, the newly developed minimum can compete 



E/B{X) = X\E2/B) + {Ei/B) + X^E^/B) + {2L^)V{X). 



(11) 
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with the one near X ~ 1. At higher density, the minimum gets shifted to X = 0. The figure in 
the small box presents the resulting values of Xmin as a function of L. There we see an explicit 
manifestation of a first order phase transition. This mechanism corresponds to reformulating 
BR-scaling|15j in a more accurate way. 

One can perform a more rigorous treatment by allowing the space dependence in XAf(^) 
(with the symmetries summarized in Table 1), which leads basically to the same physics except 
for small quantitative difference. The most essential new ingredient is that the static dilaton 
field for dense matter vanishes identically all over the space in the half-skyrmion phase. 

By incorporating the fluctuations of pions and dilaton, as summarized in Table 1, we obtain 
as before the Lagrangian for the dynamics of these particles 

L{U,x) = Lu,TT + Lu,x + I^M,iTx^ (12) 

where we have put the subscript "M" to denote that it describes the in-medium dynamics. 
Explicitly, each term can be expressed as 

LM,7r = \Gab{r)d^TTad^'TTb-\S{f)TTl 

LM,x. = P:(r-)x9'^a + ga(f)x*^ m) 

Here, Gahij)^ S{'r), V^{r), M{f), Pa{r) and Qa{^) are the effective potentials provided to the 
fluctuating fields by the background fields Um{t) and Xui^ (-^o^ the details, see Ref.jHj). 

By applying the same approximation, i.e taking the average values for the effective potentials 
over the space, we can estimate the medium effects on the properties of the fluctuating fields. 
As for the pion, we are led to a Lagrangian similar to ©, where and a take the average 
values with the additional factors {XmI fx)"^ ^"^^ (Xm//x)^' respectively. Since Xm vanishes in 
the chiral symmetry restored phase, the ratio f^/ f-n now vanishes. However, even in the chirally 
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Figure 4: In medium properties of the pion and dilaton fluctuations as a function of the baryon 
number density. In the inset figure, the pion and dilaton masses are presented in the actual 
energy scale. The numerical results on the average values of Xm ^ °f background 
matter are presented in the small figure on the left hand side. 



restored phase, with (cr) = {Xm) ~ 0' (-^) does not vanish due to the background Um{^ 
couplings to X- In Figure 4 we show the ratios of the in-medium parameters relative to their 
free-space values. 

The phenomenon discussed above is closely related to "Brown- Rho" scaling ^S]- In the 
description of Ref. ^S] , the density dependence comes solely from the change in the mean field 
X* where the corresponding change in the skyrmion structure has been ignored. Our present 
result corrects this fact and gives a precise meaning to the scaling relation of Ref. jl5j . 

We can treat the background interactions more systematically following the perturbative 
scheme we have developed. We decompose the Lagrangian into an unperturbed part, Lq, and 
an interaction part, Lj. This leads to a Hamiltonian split also in two parts 

H = Ho + Hi. (13) 

The free propagators are defined by Hq and the interaction potentials appearing due to Hj are 
summarized in Figure 5. In Figure 5, G"^(£), for example, is the Fourier transform of the local 
potential C^r): 

G"^(f) = / d^re'^-'^G^^f), (14) 

''box "'box 

where the integration is over a unit box of the crystal and Vbox is its volume. Due to the periodic 
structure of the crystal only discrete values of the momentum are allowed. 

We show in Figure 6 the diagrams used to evaluate the self-energy. Only the diagrams for 
S^^TTj, appear. The symmetry structure of skyrmion matter allows a non vanishing self-energy 
only for a = b. To first order, T,^^^ is nothing but Hi{i = 0). Since H.j^^{6) = 0, no mixing 
between the fluctuating pions and the fluctuating scalar occurs. Thus, the pion propagator for 
TTa can be expressed as 

1 = 1 (15) 

pI-p'-^^'Kpo,?) G'-(0)(p2_p-2)- 
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Figure 5: Free propagators and interactions for the pion and the scalar fields in the presence 
of background skyrmion matter. The energy-momentum conservation delta functions are not 
shown. 
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Figure 6: Diagrams used to evaluate the sclf-cncrgy of the tTq propagation up to second order 
in the interaction. Here, b runs over 1, 2, 3 and the intermediate states run over all ^ / 0. 
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where we have used that the self energy to this order is Ti^}na(.Po,P) = —p^iGaai^) — !)• The 
superscript "(1)" means that the quantities are evaluated to first order. 

Since Lorentz symmetry is broken by the medium, the general form of the in-medium prop- 
agator can be written as 



with Zf^ = {ft,s/ fn)^ and the "pion velocity" in medium is given as = Zt/Zs- Comparing 



with this, we obtain ft = fs = f iry Gaa{0)- Since Gaa{0) is nothing but the average of G'*"(r) 
over the space, our calculation thus far reproduces the naive approximations discussed above. 
To the same order, the self-energy of the scalar field is ^xx ~ -^(0) ~ ''^x- Since it is constant, 

this self-energy modifies just the scalar mass from the free value to \J M{0). 

Now, the second order diagrams shown in Figure 6 can be calculated similarly. Again, in 
spite of the tt — x coupling term in the interaction Lagrangian Ti^na vanishes so that the 
pion propagator and the scalar propagator can be simply written as 

1 I (17) 

Po-P^ - ^ntif' pl-i? -ml- 

respectively. 

The po and p dependence of T,^'^^ is not so simple as that of S*-^). By assuming that energy 
and momentum are small, we may expand the self-energy '^nant i^i powers of po and p. Then, 
we can express the propagator in the form of but Zg gets some corrections terms from the 
second order diagram. The corrections are negative definite so that the pion velocity becomes 

< 1. 

The results of our calculation are shown in Figure 7. We show the decay constants ft and 
fs in units of fj^ as a function of the density measured in a dimensional units. The inset figure 
is their ratio, i.e. the in- medium pion velocity Vt^. To second order the contribution to fs/f-K 
turns out to be small, and thus the pion velocity stays ~ 1- The lowest value is about 
~ 0.9. Furthermore, for the pions at higher matter densities, the internal propagator provides 
an extra suppression because L scales as p^^^. Once we pass the density at which the pion 
velocity has its minimum, the pion velocity increases with density and approaches 1. When 
the background matter is in the half-skyrmion phase, Xo(^) vanishes identically and so do all 
the local potentials. Thus, both ft and fs vanish. Moreover, their difference is proportional 
to the square of the potentials and vanishes faster. This fact explains the asymptotic behavior 
f 7r — > 1 in Figure 7. This result is very similar to that of Ref. 17^ found in heat bath, where the 
pion velocity approaches 1 while both the spatial and temporal pion decay constants vanish at 

T = Tr. 



4 Vector mesons 



At higher density and/or temperature we need more degrees of freedom such as the p and lo vector 
mesons. The same procedure can be repeated with the vector mesons using the Lagrangian 
The lowest energy FCC skyrmion crystal configuration can be found by requiring the symmetries 
given in Table 1 to the pion, dilaton, rho and omega fields. The numerical results are presented 
in Figure 8. 

In the irpx model, as the density of the system increases (L decreases), E/B changes slightly. 
Its value is close to the energy of a single skyrmion up to densities larger than po (L ~ 1.43). 
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Figure 7: ft and fg as a function of skyrmion matter density in arbitrary units for the chiral 
limit, m-TT = and = 720 MeV. The in medium pion velocity appears as function of matter 
density in the small box. The arrow indicates the point at which the observed ratio in pionic 
nuclei (at p « 0.6po) is located. 

because the size of the skyrmion is very small and the skyrmions in the lattice only interact at 
very high densities when their tails overlap. 

In the absence of the uj the dilaton field plays as before an important role. Skyrmion matter 
undergoes an abrupt phase transition at the density at which the expectation value of the dilaton 
field vanishes (x) = 0. 

In the TTpuJX model, the situation changes dramatically, because the to provides not only a 
strong repulsion between the skyrmions but also an intermediate range attraction. In both the 
vrpw and the npi^x models, at high density, the interaction reduces E/B to 85% of the B = 1 
skyrmion mass. This value should be compared with 94% in the irp model. In the vrpx-model, 
E/B goes down to 74% of the B = 1 skyrmion mass, but in this case it is due to the dramatic 
behavior of the dilaton field. 

In the npcox model the role of the dilaton field is suppressed. It provides only a small attrac- 
tion at intermediate densities. Moreover, the phase transition towards its vanishing expectation 
value, (x) = 0, does not take place. Instead, its value grows at high density! 

The reason for this can be found in the role played by omega in Q. In the static configura- 
tion, omega produces a potential, whose source is the baryon number density, which mediates the 
self-interaction energy of the baryon number distribution. Thus, unless it is screened properly 
by the omega mass, the periodic source filling infinite space will lead to an infinite self-energy. 
To reduce the energy of the system, the effective uj mass must grow at high density, for which 
X must grow too. Note the factor {x/ fx)"^ in the omega mass term in Lagrangian 
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Figure 8: E/B and the average values of the fields over the space as a function of L. 



5 Summary 

We have developed a unified approach to dense matter in the Skyrme philosophy, where systems 
of baryons and mesons can be described by a single Lagrangian. In our approach dense baryonic 
matter is approximated by skyrmion matter in the lowest energy configuration for a given baryon 
number density. By incorporating in it fluctuating mesons we can get some insights on meson 
dynamics in a dense medium. Our approach enables us to study this dynamics beyond the first 
order in the baryon number density. One can continue to work in this direction by incorporating 
more degrees of freedom, by improving the way of treating matter beyond the crystal solution, 
and so on. 

However, before closing the presentation, we must clearly lay down the scope of our work. 
We do not claim that the results obtained at present describe reality. The most fundamental 
problem we phase is that our "ground state" for matter is a crystal not a Fermi liquid. Our aim 
has been to assume a state for matter, given by a classical solution of a theory considered to 
be valid at large Nc, and have studied the implications for its excitations. Our work should be 
taken as representing the first step towards a more realistic treatment of a dense matter theory. 
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Table 1: Summary of the properties of the B = vacuum solution,the hedgehog Ansatz for 
the B = 1 skyrmion, the symmetries of the FCC single skyrmion crystal for pions, dilaton and 
vector mesons. 





pion & rho 


dilaton & omega 


vacuum 


U 


= 1, p2 = 


X = /x> '^M = 


hedgehog 


U = 


exp(ir • rF{r)) 
eaipfpG{r)/gr 


X = X{r) 
ojQ = uj{r) 


boundary 


m- 


= TT F{oo) = 


w'(0) = w(oo) = 


condition 


G(0) = 


= -2 G(oo) = 


X'(0) = x(oo) = 


reflection '^^^ 
(i/«-plane) 






X, ^ 


3-fold axis 
rotation''^-' 






X, ^ 


4-fold axis 






X, ^ 


rotation^^' 






FCC 






X, ^ 


translation''*'' 






Ansatz for 


u = 




X = Xm + X 


fluctuations 




= PI'''+Pl 




*U = <t>J',+ir-r, 








{x,y,z) (1) {-X 


,y,z), (2) {y,z,x), (3) {x,z 


,-y), (4) (x + L,t/ + L,2;) 
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